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> I Abstract 

in . 

^ I We compute the inclusive decay rates b -^ D~{D*~)c including lowest order QCD 

\Q • corrections on the quark legs, and compare with existing data. Unlike the short 

O . distance QCD corrections, that are of higher order, these corrections are of order 

Us- In the on-shell renormalization scheme and for ^^(mb) = 0.2 we find a cor- 
rection of —10 % to the inclusive rate computed using factorization. This gives a 
p • total rate BR{h -^ D^{D*~)c) = 8 % consistent within la with the measured value 

jL . BR{B — >■ Dfx) = (10.0 ± 2.5) %. The general formulae given here include the case 

D I of vanishing mass for the final quark b -^ D~{D*~)u. The radiative correction to 

"^ I this rate is —17 %. We show in another place that this process can be useful for the 

J^ • measurement of the CKM matrix element Vub- We also give the renormalized vertex 

^ . at the interesting values g^ = and g^ = q'^ax^ ^^^ compare with existing literature. 



0^ 



S 



LPT Orsay 99-36 
DAPNIA/SPP 99-19 
May 1999 



Unite Mixte de Recherche CNRS - UMR 8627 



1 Introduction 

The inclusive rate oi a B meson decaying into a Dj meson is obtained using the 
spectator quark model shown in Fig. 1 : 

T{B^D;X) = T{b^D;c)+Tib^D:-c) . (1) 

A main point in writing the approximation (1) is that, as we point out in ref. 0, 
the excited states D** do not lead to Ds since their dominant decays are D** -^ 
D^*^K. Moreover, the quark c dominates the inclusive rate, the quark u being CKM 
suppressed. Of course, there are also other mechanisms for producing the Dg of the 
right sign, namely through annihilation and exchange diagrams. However, as we 
discuss in detail in ref. [|1|, these mechanisms are suppressed. In this paper we will 
concentrate on the mechanism of Fig. 1, and compute the 0{as) radiative corrections 
to it. 




Figure 1: Spectator diagram for the decay b -^ Dfq. 

From the effective weak Hamiltonian, using factorization and the spectator quark 
model, a straightforward calculation gives the rates 
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T^'Hb^D;q) = ^ \V,,V:f al f^ ml^\{l,r\i) [{I - r^f - i{l + r')\ (2) 



(3) 
where Oi is the QCD short distance factor 

„ Ci c+ + c_ c_L — c_ ,^, 

Empirically, from exclusive decays, one finds ai consistent with |ai| = 1. For further 
convenience we have adopted the notation 

r = — i = ^2 (5) 

mi, rrif^ 

where rriq is the final quark mass rriq = nic or m„, g^ = rrij-,^, or m|,*, and A(a, b, c) = 

a"^ + b^ + c^ - 2ab - 2bc - 2ca. 

To compute the order of magnitude of the expected branching fractions, let us 

tentatively use the following numerical values for the quark masses (see below for a 

detailed discussion of these parameters) 

rUb = 5.0 GeV , m^ = 1.6 GeV (6) 

and rriu = and the decay constants [0 

fos = 230 MeV /d* = 280 MeV . (7) 

Using tb = 1.6 ps and \Vcb\ = 0.04, one obtains 

BR^^\b -^ D;c) ^ 3.2 % Si?(°)(6 ^ Dl^c) ^ 6.8 % . (8) 

where the superindex (0) means using the factorization results (0), (0). Within the 
assumption ([l|), this yields 

BR{B -^ D;X) = BR^^\b -> D-c) + BR^^\b -^ Die) ^ 10 % (9) 



and, for completeness, with '-yy\ = 0.08, let us give the b —^ u branching ratios : 



BR^°\b -^ D-u) ^ 2.6 X 10"^ BR^^\h -^ D*-u) ^ 4.9 x 10' 



(10) 



The naive prediction (|^) is consistent with the measured value 



BR{B -^ Dfx) = (10.0 ± 2.5) 



:iii 



The aim of this paper is to compute the lowest order QCD corrections to the process 
depicted in Fig. 1, i.e. to investigate how stable is the naive result (|^) relatively to 
QCD radiative corrections. It is important to emphasize that the calculation of 
the radiative corrections will ask in particular for a detailed discussion of the quark 
masses, that we perform in Section 6. 
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Figure 2: Diagrams leading to the production of Ds mesons. 



Let us first notice that the leading logarithm corrections cancel at order a^. Indeed, 



"4m) 



where d^ = — ^, 



from the well known expressions at leading order c±(/i) — ^ (m ) 
(i_ = ^, one finds for the expression (^) ai = 1 + O^a"^). The correction is at most 
of second order, and this explains why the combination of short distance coefficients 
ai is predicted by the theory to be close to one. This agrees with the measurements 
on class I exclusive decays 0], that yield an empirical |ai| = 1.00 ± 0.06 when one 
uses SVZ factorization to calculate the matrix elements. 

Figure 2 depicts the lower order QCD corrections. The contribution of Fig. 2a 
corresponds to gluons going from the lower to the upper vertex, they are UV divergent 
and contribute to the anomalous dimension. Figs. 2b, c are UV convergent if one adds 
the corresponding self-energy diagrams on the quark legs : this corresponds to the 
well-known fact that currents have no anomalous dimension. The fact that, as we have 
seen, the contribution of Fig. 2a vanishes at 0{as) follows from the color neutrality of 
the DI*\ Indeed, the W produces a gg in a color singlet. If the latter emits only one 
gluon which goes to the other vertex, and does not reabsorb any gluons, it becomes 
a color octet and cannot end into a D'f'. Of course, this vanishing at 0(0;^) holds in 
general, not only for the leading a^ log ( ^^ ) term in the combination ai (^. 

Since 0{as) leading log QCD corrections vanish, we will have to consider 0{as) 
corrections not containing the coefficient log ( ^^ | , i.e. we will have to consider also 
the UV finite radiative corrections affecting separately the lower and upper quark legs. 
These corrections are IR divergent, and have to be combined with the Bremsstrahlung 
graphs on the same leg (Figs. 2d,e) to obtain the IR convergent physical result. 
At the upper vertex, however, all the gluon corrections of the type 2c, e (and the 
corresponding self-energies) have to be included in the vacuum to D^*^ matrix element 
represented by the decay constants /d^ or fu*. Therefore, if we use a consistent 



evaluation of the latter, we will not have to bother about gluons at the upper vertex. 

In the following, we present an estimation of the 0{as) radiative corrections at 
the lower vertex {b-q quark line), for arbitrary q quark mass and g^ carried by the 
current. Of course, the interesting cases are mq = rric or tUu and q^ = m\,^ or rn^)*. 
Calculations of radiative corrections at the order aimed in this paper have been made 
in the past by Guberina, Peccei and Riickl 0, that computed the 0{as) corrections 
to the total rate h -^ q'iq'2^3 in the case of vanishing masses for the final quarks. The 
result is the sum of two contributions : the QCD translation of the Berman, Kinoshita 
and Sirlin calculation of QED corrections to // decay and the 0{as) correction to the 
ratio R in e+e~ annihilation. Guberina et al. used, as we do in this paper, naive 
dimensional regularization (NDR) |^ for both UV and IR divergences. From their 
intermediate results for what is now for us the lower current in Fig. 2, we can compute, 
saturating with the tensors /|, q^q^ (for the D^) or fi)*{—g'^'^q'^ + q'^q") (for the sum 
over the polarizations of the D*), with g^ = ?Ti|, or nij-,*, the total rates r(6 -^ Dju) 
and r(6 -^ D*u), since the u quark mass is approximately massless. We have done 
it, but we need to make the general calculation for a final massive quark q. In the 
limit nig — i> we recover the result that we have obtained from the intermediate 
stages worked out by Guberina et al. Here we will expose and discuss the general 
calculation. 

Radiative corrections for weak processes taking into account the unequal quark 
masses have deserved much attention in the last years. In the nonleptonic three 
quark decay b -^ qiq2q3 with arbitrary quark masses, the calculation of the radiative 
corrections was performed by Hokim and Pham H]. At lowest order, the current 



vertex had been computed by Gavela et al. and Halprin et al. ITH], and by Paschalis 



and Gounaris and by Schilcher et al. |T^. The radiative correction to the semileptonic 
heavy quark decays was performed by Cabibbo and Maiani and by Nir |TB|. Recently, 
with the aim of applying the results to Heavy Quark Effective Theory, the vertex has 



been reconsidered by Voloshin and Shifman |T^ and by Neubert |]T5 . 

The paper is planned as follows. Although these are, as the folklore says, "straight- 
forward but tedious calculations" , we will give some details, helped with Appendices 
in order to make the discussion simpler. In Section 2 we summarize the calculation 
of the renormalized vertex and in Appendix III we compute the limits at g^ = and 
at g^g^x to compare with the existing literature. In Section 3 we give the corrected 
two-body rate. In Section 4 we summarize the calculation of the Bremsstrahlung rate. 
In Section 5 we discuss the analytical results and the necessary cancellations in -^^^ 
and of the mass singularities, and we deduce the rUg —* and q^ ^ limits of the 
final finite correction. In Section 6 we give numerical results for several interesting 
cases, along with a discussion of the values of the quark masses. 

2 Renormalized vertex 

We use NDR regularization and follow the Feynman rules and notations from the 



review article by Aoki et al. on the Standard Electroweak Model |T^. In particular, 
Dirac algebra is performed in D dimensions in Minkowski space and we work in the 
Feynman gauge. 

A long but straightforward calculation gives, for the renormalized vertex 

7,(l-75)AL(r,0 + — (1 + 75) 5+(r,0 + — (1 + 75) C+(r,0 



mh 771), 



(12) 



where the coefficients of the different Dirac structures are 



Mr,0 = l 


- 


D-3 {D- 
D-A ' 2{D 


l)(l + r^-4) 
-4)p-3) 


r2 ^ l-r2 1 




B^{r,0 = lh 




CAr,0 = -l 








Anir.O ro-,'^ 




B-{r,0 = - 


1 
2 


r 


[ih'- 


D 
D 


J/3 + /.] 





1 - C + r^ 



C_(r,0 = 2 ^(^2-2/3 + 14) 



(13) 



In the second term of the expression of Ai{r, C,) we recognize the vertex counterterm, 
that we deduce from the Ward identity in Appendix I, and agrees with the resuh 
of [^. The quantities Ii{C,,r) are the integrals 
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D/2~3 
D/2-3 



(14) 



whose expressions we give in Appendix II. 

As a partial check of our calculation, we have computed the renormalized vertex at 
g^ = and at q^ = q'^^x- We give our results in Appendix III. The renormalized vertex 



for arbitrary masses and q had also been computed by Paschalis and Gounaris [12 



using as infrared regulator a gluon mass A, instead of dimensional regularization as 



we use here. At g^ = it had been computed by Gavela et al. |T3| and by Hokim and 
Pham 1^. For the IR finite form factors we agree with their resuh. At q^ = q^^xi ^^ 
agree with Paschahs and Gounaris [O, Voloshin and Shifman |T3], and Neubert [113 



3 Two-body decay rate at order a^ 



After a rather long calculation, one finds, for the two-body decay rate at 0{as) '■ 

4 a, 

e 



-pTwo—body p{0) 



2Al + 






B+ + 



[l-i-r 



2\2 



+ 



4r^ ^(l + ^_^2^2 



C^ 



(l-r2)2-^(l + r2) 
r[(l-r2)2-^2] 



Ar + 



B^ 



:i-r2)2-^(l + r2) 

-nil 

(l-r2)2-^(l + r2) 
where Al, • • • C_ are understood to be functions of r and ^, and the function g{D) is 



(15) 



r 



6-D 



VTTTl 



given by 

9{D) 
Expanding in powers of -D — 4 one gets finally 



r(l) 



2\ D-A 



(27r)«-4ff(47r)^-4 r (^^] \ yU^ 



A(l,r^O 



(D-4)/2 



(16) 



^Two-body _ p(0) 



£>. 



D, 



1 + 



4 a. 



3 TT (1 — r 



2^2 



2 (Ai + Di) + 

4r^ 
^(l-r2)2-^(l + r2) 



:i-r 



2-\2 



e^ 



;i-r2)2-^(l + r2) 



-e(l + r2 

5+ + 



2n2 



;i-e-r2) 



A 



R 



r(l + .^ — r 



;i-r2)2-^(l + r2) 

2\2 



c. 



(l-r2)2-^(l + r2) 



fi_ 



r[(l-r2)2-^2] 
;i-r2)2-^(l + r2) 



C_ 



(17) 



where DL{r,C,) comes from the expansion of g{D) and the terms of order -j^^ in 



^L(r,0 



D, 



l--{l-i + r^)h 



-2 log(27r) + 7 - 1 + log 



Tim 



IX' 



"^ f-log[A(l,r2,0 



'1^ 



In this expression I2 means the hmit I2 ^ that can be read from the formula in 
Appendix II. Analogously, the other terms become 



A.^-J 



1 r 1 



2 ^D-4: 

+ (1-^ + 20/2} 



3 - Ji - (1 - e + r')/2j -I1-2 + - log(r) + (1 - r^)h 



+ 2 



C4 



-2(^3 + /4) 



S_ = -^ r (2/2 - 3/3 + /4) 



(19) 



where the limit D ^ A has been taken wherever it is possible, i.e. everywhere except 
in the first term of Al. 

4 Bremsstrahlung rate 



From the expression of the real gluon emission amplitude 



^7^(1-75) 



Ut^ ' ^ ; -^TTT 7m (1 - 75) 



57^(1-75)0 



(A) ^^ 

2 



2 m2 - (p„ + A;)2 
one gets the color and spin averaged rate for Dg emission : 



(20) 



I A/<Bre»n|2 



1 ^2 2 f-2 fi V 
o *-^ Qs JD Pd Pd 
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u-f^, (1 - 75) 2 / r\^ 7a& + W7a o / —r\^ 7^ (1 - 75) b 



1^1 - {Pb - kf 

- (^ ^ h-f^ + ^b 

^7^(1-75) 



7"6 + ui 



^l - ipu + ky 



lu{^ -l5)b 



.(21) 



ml - {pb- ky ' " ' "" ' ml- {pu + ky 
Performing the Dirac algebra and putting quarks and gluons on-shell, one arrives at 
the expression for the rate : 
2 1 



^Brem 



3 (27r; 



2D-3 



G' gl fn ml 



[l-r 



2\2 



e(l + r^ 



j-2,0 



1 



mt 



[l-r 



2\2 



ai + r^ 



-1,0 



+^4(l + rV'° + ^2 



ml 



mt 



l-r^f^^{l + r^ 



rO,-l 



+ r^ 



'1 — r 



2\2 



e(l + r^ 



r0,-2 



_i_2(l + r2)/-i'i-i-2(l + r2)/i'-i ^^— ^ ^^, .2^2 



ml 
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l+r'-O l-r^)'-e(l + r' 



-1,-1 



(22) 



where J'"'" are the phase-space integrals in D dimensions : 



j^,n ^ ^i^+n)(D-4) / fZ^ fZ^ fZ]Jl ^^(p, _ p^ _ p„ _ A;) (p, ■ A;)'" (p„ ■ kf 



2pl 2p% 2A;0 



that are functions of r and ^ and whose expressions are given in Appendix IV. 



5 Discussion of the analytical results 



5.1 -D, rate 



(23) 



It is convenient to reorganize the expressions of the two-body and Bremsstrahlung 
decay rates in order to check the necessary cancellations, namely the -^^ poles and 
the mass singularities. We can reorganize the expression for the two-body rate in the 
form : 



-pTwo—body F^'^-' J 1 _|_ 



4 a, 

3 TT 



X + Yl + 



{l-r^f-( 



;i-r2)2-^(l + r2 



B+ 



11 



+ 



;i - e - r 



2\2 



'1 — r 



2A2 



ai+r^ 



C+ + 



4^r 



1 — r 



2^2 



e(i + r2 



A 



i? 



, r(l+e-r^)^ ^ , r[(l-r^)^-e^] ^- 



;i-r2)2-^(l + r2 



;i-r2)2-^(l + r2) 



(24) 



where 



X 



l-^ + r2 , /l+r2-e + v'A(l,r2,0 
log 



VA(l,r2,0 Vl + r2-e-v/MM^. 



X 



2 log(27r) + 7 + log 



Trr/it 



(25) 



and the new quantity Y^ is finite for D ^ 4, given in Appendix V. Analogously, we 
can rewrite the Bremsstrahlung rate in the form : 



-\Brem 



<0) 



4 a. 



— i- n - 

3 n 



D. 



2K~''" + 2K 



0,-1 



X + ir-2'O + r^ir"'-^ _ (^1 + r2 - ^) J^-i'-i 

4(1 + r2) 



K 



0,0 



(l-r2)2-^(l + r2) 
2(1 + r2) ^_,, 2(1 + r2^ 



ir 



-1,1 



.i^i.-i 



(26) 



(l-r2)2-^(l+r2) (1 _ ^2)2 _ ^(1 + ^2)- 

where X is the quantity defined by (|25| ) and the expressions /i'"*''^ are given in Ap- 
pendix V. We observe that X, that contains the j^ terms and the terms in log ( ^^) , 
cancels between the two-body and the Bremsstrahlung rates, as expected. To check 
the cancellation of the rest of the mass singularities, we need to perform an expan- 
sion in powers of rrig or in powers of r. From the expression (^) and Appendix II we 
obtain, for r — > : 



-pTwo—body p(0) j -i , 



Ds 



4 a, 

3 TT 



D(r,0 



r^O 



+ 



Fv{r,0 



r-^O. 



(27) 



[D{r,^)]r-^Q contains the ^^^ and the divergent terms as m^ — > 



+ log(r) 



r^O 



2 - 2 log 



'l~C 



D-4 



- 2 log(27r) + 7 + log 



1^. 



-- + 21og(l-0 



+ log^(r) 



(28) 
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and [Fy(r, ,^)]j,_>o is the surviving finite piece : 

[Fv{r,OUo = -3 - y - 31og2(l -0 +41og(l -0 + ^ log(l - 
+^p(l-0 + logeiog(l-0 • 



(29) 



The Spence function Sp{z) is defined in Appendix II. Analogously, we obtain 



D. 



3 71 



(30) 



[FB{r, oUo = -f + y + y ^°^^^ " ^^ + ^°^^^^ ^°^^^ ~ ^^ + ^^^^^ 

-31og2(l-0 + Y^log(0 • 



(31) 



We observe also that the singular terms in log(r) and log (r) contained in [D{r, 0]r^o 
cancel among the two-body and Bremsstrahlung rates. 
Then, it follows, for rriq -^ 0, the total rate : 



Two— body j^ -pBrem p(0) 



Ds 



+ r 



Ds 



Ds 



1 + 



4 a. 



3 71 



9 tt' 



2SpiO - log(e) log(l - 



+ ^ log(l -O-l log(l -0-^3^ log(0 



(32) 



that gives, at g^ = 0, i.e. ,^ = : 



Two-body . -pBrem _ 1^(0) 



Ds 



+ r 



Ds 



Ds 



■ 4 a, /5 _ TT^^ 

3 TT U 3 , 



(33) 



5.2 D* rate 



The calculation of the D* rate follows along the same lines, with the replacement 



f2 _2 „!yi^ 



Jd Pd Pd ^ Jd'^d* Z^^ ^ — Jd*Pd*Pd* Jd'^^d*9 

X 



(34) 
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The difference between the D* and Dg rates is thus a term proportional to m|)*, i.e. 
q^ or ^. In the hmit .^ ^ one must recover the same decay rate. We obtain, for the 
two-body rate : 



Y^Two—body 
^ D' 



+ 



A(l,r^O 



rA(l,r^O 
(l-r2)2 + ^(l + r2-20 

and for the Bremsstrahlung : 






A(l,r^O 



I2ri 



B, 



B_ 



l-r2)2 + ^(l + r2-20 

rA(l,r2,0 
;i-r2)2 + ^(l + r2-20 






(35) 
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.2 7^0,-2 
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+2ir 



■0,-1 



4(1 + 



:i-r2)2 + ^(l + r2-20 



is: 



0,0 



2(l + r2 + 20 



(l-r2)2 + ^(l + r2-20 



is: 



-1.1 



2(l + r2 + 20 



K 



1,-1 



(36) 



(l-r2)2 + ^(l + r2-20 
Again, as expected, X, that contains the -jj^ terms, cancels between the two-body 
and the Bremsstrahlung rates. As for the mass singularities, we find in the case of 
the Dl exactly the same expression for [D]r^o that cancels among the two-body and 
the Bremsstrahlung rates. The finite result is, in this case : 



sTwo—body . -r^Brem 



-pi u 



Dt 
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Dt 



4 tts 

1 + 7T- 

3 n 



71 



2---2Sp(0-log(Olog(l-0 



2(1 + 20 
that gives, at .^ = : 



log(l - 



e(l-e-2e^) 

4(1-0(1 + 20 (1-0^(1 + 20 



log(0 



(37) 
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Dt 



■ 4 a, /5 _ 7r2 ^ 

3 TT U 3 , 



(3^ 



i.e., the same correction than for the Dg, as expected. 
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6 Numerical results 

To give the numerical results, let us parametrize the rate, including the QCD 
corrections, under the form : 



r (6 ^ D;q) = rw {b ^ D;q 

and analogously for the D* 



6 TV 



(39) 



r (6 ^ Dl-q) = rW {b -. DTq) \l + \-ri* (iDt.r) 



(40) 



3 TX 

The functions ri{S,i r) and ri*{^, r) can be read respectively from the finite sums of the 
two-body and Bremsstrahlung corrections, formulae (p^, ( pG]) and (|35D, (pGl). The 



particular values ^(■Cds)'") and ri*{C,D*,r) correspond respectively to both functions 
for g^ = m|)^ and mj^*. At fixed r these are slowly varying functions of C,- For r = 0, 
these functions are given by the relatively simple expressions ( |32D and (|37D, and for 
r = Tf, = ^^ we use the full expressions. The functions decrease monotonously with 
^, the dependence on g^ is mild, and the values at the interesting values of low q^ (i.e. 
^ < 0.2) do not differ significantly from the g^ = value. For ^ -^ imax there is a 
logarithmic singularity, of the form log(l — ^) for r = 0, smoothed out by the phase 
space V^^^{^rnax) = 0. The functions ?7(^,0), 'r]{^,rc), i]*(^,0) and i]*(^,rc) decrease 
by about 50 % between ^ = and C, = 0.6 ^max- The correction depends weakly on 
the ratio r = ^^ in the expected range of quark masses. The radiative corrections 
depend on the ratios (,£), ^d* and r and on as{f^)- We take fi = mi, and we adopt the 
value as{mh) = 0.2. 

For the quark masses, we take pole masses from a fit to the semileptonic decay rate 
b -^ ci^Ui with QCD corrections at one loop [|l^], to be consistent with the same order 

15 



that we compute here. The semileptonic decay rate reads, in this approximation : 

(41) 
where the phase space fpsii^) and radiative correction fncij) functions depend on 



V{h -^ d V,) = ^^ \Var fps{r 



1 + TT - fRcir) 

6 TX 



r = ^^ and are given in |]T6|. Setting r = 0.3, we obtain, from the semileptonic 



branching ratio 11 %, 

mh = 4.85 GeV m^ = 1.45 GeV . (42) 

The mass difference rrih — rric = 3.40 GeV compares well with the value m^ — rric = 
(3.43±0.04) GeV obtained in the 1/m expansion of the Heavy Quark Effective Theory 



1% . Therefore the pole masses that we choose (^2D, at one loop, seem reasonable. 



With these values, and m„ = 0, the results are the following, for the u-quark : 

4 ot 

--v{^Ds,0) = -0.168 

6 TT 



and for the c-quark 



^-V*UDt,0) = -0.159 (43) 

3 TT \ " / 



4 a 

-—vUDs,rc) = -0.096 

6 TT 



^-^V*{iDt.rc) = -0.108 . (44) 

Then, using these values, the corrected branching ratios will be, from (^ and (|T0|) : 

BR (h -^ D-c) = 2.6 % 

BR (b -^ Dl-c) ^ 5.4 % . (45) 



Our conclusion is that the sum including radiative corrections 

BR{h^D-c)+BR{b^ D*-c)^8% (46) 
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is still in agreement within la with the measurement (|TT]). Also, the moderate radia- 
tive correction to the processes b -^ D^*^~u reinforces our argument about using the 
spectrum B —>■ DjX^ in the measurement of the CKM matrix element Vub ^■ 

A number of remarks is in order here. First, from the relation at one loop between 



the pole mass and the running MS mass 



m = 771(171) 



4 ajm) 



3 71 

at first order in as, from (^I]) one obtains : 



(47) 



1927r3 1^^"' fps{r)<l + - ^ 



Tib -. ciu,) = " 'ZT' IK.P fPsir) 1 + 1 ^^^^ [5 + fncir)] (48) 



and in formulas ( p9D and (|iOD when [mb{mi,)]^ is substituted to ml, a term |^ x 3 
has to be added. Taking as{m{,) = 0.22 we get from ( ^7D 

rhbirhb) = 4.43 GeV . (49) 

The values of the pole mass (|4^ ) and of the MS running mass (^91) are respectively 
smaller and larger than the values recently quoted in the literature from the analysis 
of semileptonic b decay, the reason being that a partial resummation of higher oder 
diagrams is made that enhances the radiative corrections by roughly a factor 2 (see 



for example [T^). The interest of considering the MS mass is that the series is Borel 



summable, while using the pole mass, there is a renormalon ambiguity, cancelled by 



another renormalon ambiguity in the pole mass. Ball et al. |17| quote, as central val- 



ues, rrih = 5.05 GeV, rric = 1.62 GeV, and mb{mb) = 4.23 GeV, mcirric) = 1.29 GeV, 
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leading to consistent results for the semileptonic rate in both the MS and on-shell 
schemes. 

Concerning the radiative corrections at higher orders in the processes h -^ Dy^"c, 
one part of the radiative corrections has the same topology than in semileptonic 
decays (simply at a given value of q^ instead of integrating over g^), and we could 
expect that these radiative corrections would be similarly enhanced by higher orders. 
This is by the way what happens at order as '■ the correction that we obtain for 
b -^ D^*^^c is very close to the one obtained at the same order in the semileptonic 



case 



16| . This would imply, e.g. in the on-shell renormalization scheme, a larger 



radiative correction by a factor 2 |T^, but consistently also a larger rrih = 5.05 GeV, 
leading grosso modo to the same results (^) and (^61) . We must keep in mind however 
that, already at second order in a^, we could have another type of corrections, absent 
in the semileptonic decay, that break factorization (e.g. two or more gluons linking 
the Ds to the quark legs). We can only hope that these corrections will be small, like 
it is the case for the short distance QCD coefficient (H), that empirically is very close 
to 1. 

Work remains to be done, in particular the calculation of the QCD-corrected spec- 
trum b — > Djc to be compared with the measured spectrum |jT8| oi B —>■ DfX. It 



would be interesting to check in particular if the same function describing the 6-quark 



Fermi motion ||T9[ fits the semileptonic spectrum and the inclusive Dg one as well. 



Acknowledgements 

The authors acknowledge useful discussions with J. Charles, and partial support 
from the EEC-TMR Program, contract N. CT98-0169. 



18 



Appendix I. On-shell renormalization, Ward identity and ver- 
tex count erterm 

The bare self-energy of a quark of mass m writes 



with 



A{p') 
B{p') 



3 ^' 2^7r«/2 



"3 9s 2D^D/2 



A-D 



r 



2 
A-D 



m 



D-4 



D^ 



A-D D 



m 



D-4 



2D 
D^2 



m F 



2;^ + l;« 

A-D . D 



2 "'^'2'^ 



To proceed with the on-shell renormalization it is convenient to expand S(p) in powers 
oi "p — m, that gives : 



S(p) = -77 9 



1 



3 ^' 2^7r«/2 



A-D 



m 



D-i 



2{D + l) 1-D , 

m + ^; r p 



D-3 



D-3 



+ ■ 



The unknowns a, b in the self-mass counterterm 



a{p) = af + h 



will be fixed by the on-shell renormalization conditions 10 



u 



(p) S(^)(p)| =0 S(«)(p)m(p) 



5/=m 







where u{p) is a solution of the Dirac equation with mass m : {'p — m) u{p) = 0, and 
5](^)(p) = S(p) + a{p). These conditions yield 



4 2 1 j^f^-D\^D-, 



^(P) = -^9s ^o^r>,2 



m 



2{D + i) 1-D ; 

m + ^^ p 



D-3 



D-3 



On the other hand, as we use NDR, that preserves chiral symmetry, both the self- 
energy a{p) and current A^^ counterterms will satisfy separately the Ward identity 
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(an overall factor ^^V^fe is understood, where g is the weak coupling) |T0| |[TT[] : 

ip, - PtV Af = [a''(j9,)T(+) (1 - 75) - (1 + 75) T(+) a\p,)] - tMw X'' 

where A''* is the counterterm of the coupling of the unphysical Higgs and T^^^b = q. 
One gets immediately 



r'.^,j_Ur^«-^ 



3 ^^ 2«7r^/2 V 2 ; {D-3){D-A) 



m?-' + mf -^ 



7m (1 - 75) • 



Appendix II. Vertex integrals 

An expansion is made of the vertex integrals /j(^, r) up to first power of D — 4 in 
the cases in which this is necessary (Ji and I2) '■ 



h = l 



D-A 



1 



2 + 



log 



^Kl.T\^) /l + r2-e + VA(l,r2,0 
^og[r ) + — log 



2e 



l + r2-e+VA(l,r2,0 



,l + r2-e- jA(l,r2,0 



1 



A(l,r2,0 Vl + r2-e-VA(l,r2,0. 
1 



l + -p-4)loge 



+ -(/}- 4) 
2 -/A(l,r2,0 



-25j9 



2A/A(l,r2,0 



,l-r2 + e + v/A(l,r2,0 



2Sp 



2A/A(l,r2,0 



,l-r2-e + VA(l,r2,0 



■^W 



l-r2 + ^-VA(l,r2,e)\ 1 Jl-r^-i-\l\(l,T\i) 

+ :t lo; 



2 \i-r2 + e + VA(l,r2,0y ^ " ^i _ r^ - ^ + ^/A(l,r2,0, 

/l 
+ log' 



2, l-r2 + e-v'A(l,r2,0 



2e 



-log' 



2^ 



log 



+ log' 



2,l-r2-e-VA(l,r2,0 



2e 



2^ 



is = iog(r j = log = 

2e 2eVA(l,r2,0 Vl + r2-^-./A(l,r2,0, 
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1 l-r^-i, , 2N 

h=- + TTTK-^ log(r2) 



i 



2e 



:i - ^2)2 + e(e - 2) j^ / 1 + r2 - e + V^All, r\ 



2ex/A(l,r2,0 



l + r2-e-VA(l,r2,0 



The Spence function Sp{z) is defined by 

Sp{z) 



'*!25W 



t 

oo n 



In the hmit q^ — > Q'max' since S,', 



h{Ua.,r) = l-{D-A) 
1 



z| < 1 

(1 — r)2 and A(l,r2,,^„ 
r 



0, one obtains : 



1 



1 — r 



log(rj 






l + iD-A)[l 



I 



1 — r 



log(r; 



1 



1 



1 — r 



log(r) + 1 



1 — r 
hiUax, r) = j^-^ + jY^ log(r 

Another hmit used in the text is g^ or ,^ ^ : 



/i(0,r) = l 
h{0,r) 



D-A 



log(r2) — log(r2 



1 — r^ 



1 — r^ 
1 



log(r2)--(D-4)- ^log2(r2) 



+ 



r 



4 



log(r^ 



I /y -Li I I rf A \Zi 



(l-r2)3 



log(r^ 



2(1 -r 



2^2 



Appendix III. Renormalized vertex at q^ = and q^ = q^^^j. 

Vertex at g^ = g 

Defining the quark couphngs by the expansion of the renormahzed vertex (g = 

Pb-Pu=p-p') ■ 

\(m/ ' \ , 9m . u , 9t . u , 9s . 9p 



2mh 



2mh 



2mh 



2mh 
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we find, at g^ = : 



4 a, f 5 1 l + r2, „^„^ _3 _ 7^2 _ gr + 2(1 + r)^ 



3 77 1 8 "^ 16 1 - r2 



16(1 -r2) 



log(r^ 



1 

+2 



1 , '•— ' 

h lOGT 

D - 4 ^ \ /i- 



^)+7-21og(2vr) 



2 + 



1 + r^ 
l-r2 



log(r 



4 a, f 5 , 1 l + r2, ,, ^, _3 _ 7^2 + 8r + 2(1 - r^2 
16 l-r2 



^A = 1 + 77 — < -- + 
3 TT 



2/„2\ 



log (r 



1 

+2 



TTrri' 



D-A 



+ logl^l+7-21og(27r) 



2 + 



16(1 -r2; 



log(r^ 



1 — r'' 



■ log(r 



4 Os 
9m = - — 

3 TT 



1 - r r(l -r) , , ^x 
+ 777^ ^ log(r2) 



2(1 -r2) 2(l-r2)2 



9t 



4 a, 

3 77 



1 + r r(l + r) , , 2x 



2(1 -r2) 2(l-r2)2 



4 a^ 
9s = 1: — 

3 TT 



(l-r2) 



(1 + 2r - 2r2 - r^) + (3 + r - r"^ - ?,r^) \og{r'^) 



2A2 



2 1-r 



9p 



4 a, 

3 TT 



(l-r2) 



2A2 



'l-2r-2r2 + r^) - 



2(l-r2)3 



(3 — r — r + 3r ) log(r 



The finite couplings qm, 9t, 9s, 9p agree with former calculations by Halprin et al., 



and by Gavela et al. [13|, and by Hokim and Pham H]. We find infrared divergent 



results for gy, 9a, as we found in |T^ with a gluon mass A as infrared cut-off. In our 
expressions for gy and g^ the last term, that vanishes for r — >■ 1, cancels when the 
two-body decay rate is added to the Bremsstrahlung rate, just as it happened with 
the A regulator. In the equal mass limit r — i> 1, we find, keeping fixed the infrared 
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cut-off, the expected result : 

4 a. 



Af(p',p)=7^(l-75) + 



3 TT 



1 1 1 

7^1/1 75 



2 2m'^-^^^'"^") + fi^(^-~^'J^^ 



7 1 
6 2m 



Vertex at g^ = q. 



2 
max 



From ([T3|) and the hmit of the integrals (14) at q^ -^ Qmax given in Appendix II, 
we obtain 



Afb',p) = 7Jl + ^f 



1 + r 
-1 — -TTTz T log(r) 



+ 



-7m75 i 1 + 

4 as j 1 
3 TT 1 2m6 



4 a. 



3 7r 



^0-ui.g 



3 1 + r 
2 ~ 2(l-r 
1 



2 1-r 



log(r) 



?M 



2mh 



2(1 -r) 
1 + r 



log(r) 



+ 



2mf 



i(^^Luq 75 



1 1+r 

+ 777-^ TTT log(r) 



1-r 2 1-r 



1-r 2(1 -r)- 



log(r) 



2mb 



2(1 + r) 1 + r^-lOr 



1-rl 



2(1 -r) 



log(r) 



in agreement with Paschalis and Gounaris []T2| . In terms of four- velocities, this reads 



Ai/^)(P,P')=7Jl + 



4 a, 

3 IT 



1 3 1+r 



-7m75 |i + 



4 a, 

3 TT 



3 1+r , , , 
-2 — : z log(r) 



^4 a, j 

+3vr- 



4 1 -r 
r 



+^^/.75 



2(1 -r) 2(1 -r)2 

3 + r r(5 — r) , , , 

log(r) 



log(r) 



+ v' 



+ 



+ 



2(1 -r)2 2(1 -r) 



+ ^75 



2(1 -r) ■ 2(1 -r)2 
r(l + 3r r(5r — 1) 



2(1 -r)2 2(1 -r) 



log(r 
log(r) 



that is also in agreement with the calculation of Neubert (first reference [|T5|). Finally, 
using the Gordon identities for unequal masses 



uip'hf, u{p) = - u{p) 






u{p) 



u{j)' 






75 u{p) = 
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one obtains the form 



AfHp',p)=7.(l-75) 

3 3 1+r 
2 ~ 4 1-r 



4 a 

+ 3V1^- 



log(r) 



7m75 



3 1+r , , , 
-2 : log(rJ 



4 1 -r 



1 1 

2 2 



1 
2 



io^,{v- -v'-)^-ia^,{v- -v'-)^^ 



3(1 + r) 



2('l-r^ 



+ 



3r 



l-r) 



log(r) 



+ 2^^-^^ 



1 + r 



2r 



log(r) 



4( 



2 v^^ ~ ^mJ ^5 



2(1 -r) 2(1 -r)2 
3r2 + 2r + 3 4r(l + r 



2(1 -r) 



+ 



2(1 -r)3 



log(r) 



The vector and axial vector couphngs have been written down by Neubert [|T^], and 
we agree with his result. In the equal mass limit r ^ 1 we recover, as we must, the 
result of equal masses at g^ = given above. 



Appendix IV. Phase space integrals 



Two-body phase space 



2vl 2p% 



Q— S {Pb-PD- Pu) 



n 



{D-l)/2 



nib 



D-A 



u 2^-2 r (^) V /^ 

Three-body phase space integrals 



A(l,r^O 



(Z)~3)/2 



Tm,n {m+n){D— 4) / 

-/i j 



2p0 2^^ 2F~ 



^ {Pb - Pd - Pu - k) {pb ■ A;)" {pu ■ kf- 



''"-^'4{u^v/^^^^ 
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+ 



log 



Tvmt 



fi^ 



+ 7 - 2 + 3 log VA(1, r2, - log(r) - log [e^^ 



1-e + r" lo, 



l + r'-^ + J\{l,r^,0 



2r 



1 + e-r" lo, 



^/H^,r\0 






r0'-2 = TT^ 



1 1 



mt r 



+ 



log('^')+7-2 + 31og(v'A(l,r2,0 



-log(r)-log(e^/2 



\/A(l,r2,0 



2(1-0 log 



l + r^-e + VA(l,r2,0 
2r 



1 - e - r ) log ( ^^ 



' ' = TT 



1 



m^ I D - 4 



log 



l+r2-e + VA(l,r2,0 



2r 



+ 



Trmt 



21og 2 

/l + 

+ 3 log 



+ 27-2 + 4 log(r) + log 



^+r^-e+yA(l,r2,0 
2r 



l-r^ + e + VA(l,r^O 
2e 



log 



l+r2-e+VA(l,r2,0 



2r 



+ log(r) log I ^::^^ ' ~ T ^ 



2e 



l+r^-e-VA(l,r2,0 
2r 



-Sp 



l-r' + ^-J\il,r\0 



Sp 



l+r^-e-VA(l,r2,0 



/-''° = vr2<|^V'A(l,r2,0-^'log 



l+r'-^+J\{l,r^O 



2r 
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-^log 



'l-r' + ^+^\{l,r^,0 



2^1/2 



^"' =^ -o\/Ml,r^o + (l-Olo? 



l + r^-^ + J\{l,r^,0 



2t 



-^log 



2^1/2 



1 ri 



4 U 



m^- - l-3r^ + 5eh/A(l'^''0 



+.Mog(ii:!^i±V^i^^l_e(2-2.^+aiog^^-^^^^^^/^^^ 



2r 



2^1/2 



/^'- =^'<\ {-\ (3 - r^ - 5^) yX(i:^ 



+2 



;i-0' + 2r2e log 



l + r2-e + ^/A(l,r^0 



2r 



-2^ 2-2r2-e log 



2^1/2 



ro.o = vr^ m^ 



^ {^(l + ^' + e)^^!,^^^ 



-2r2(l-01og 



2r 



- 2^(1 - r^) log 



2^1/2 
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Appendix V 



Here we define necessary expressions that enter into formulae (P^, ( PBD and (|53|), 
(B3), that are not given in the course of the text. For the two-body decay rate we 
have the new expression : 

1 - ^ + r^ 



Yr 



log 



l + r^-i+\Kiy.i) 



V/A(l,r2,0 Vl + r2-e-.yA(l,r2,0 



-l + -log[A(l,r^O 



1 - ^^ - ^ , . 9x 

-2 + ^^ log(r2) 



-^ log 



,l + r2-e-VA(l,r2,0 



+ -log(r) . < -log Clog , 



—Sp 



2jX{l,r^0 



--log 



1 , 

+ 2I0I 



l-r2 + e + v/A(l,r2,0. 
2, l-r2 + C-^/A(l,r2,0 



+ 5p 



2VA(l,r2,0 






1, 
7 log 



2,l-r2 + e-v/A(l,r2,0 



2e 



2 °^ I 2e 



— log' 
2 ^ 



l-r2 + e+^/A(l,r^0 



+ 2lo^ 



, 1 - e + 2r% 
+ , =log 



2^ 



2^ 



A(l,r2,0 Vl + r2-e-VA(l,r2,0 
log(r ) H ;_ log 



2e 



A(l,r2,0 Vl+r2-e-v'A(l,r2,0 



And for the Bremsstrahlung : 



K 



-2,0 



2 + 3 log VA(1, r2, e) - log(r) - log(e'/2) 



log ;7-^^ ; log 



A(l,r2,0 



2r 



VA(l,r2,0 



2^1/2 



27 



K' 



0,-2 



2(1-0 
VA(l,r2,0 



-2 + 3 log JA(1, r2, - log(r) - log(e'/2 



log 



l+r2-e + V^(l,r2,0\ 1-e-r^ , /l-r2 + e + x/A(l,r2,0 



2r 



V/A(l,r2,0 



log 



2^1/2 



K 



-1,-1 



1 , /l+r2-e + JA(l,r2,0 

log I :-^ I X 



V/A(l,r2,0 



2r 



-2 + 41og(r) + log I ::^ | + Slog 



2^ 



2r 



+ 



1 



^K^r^O 



log(r) loe 



l-r2 + e + xA(l,r2,0\ TT' 



2^ 



-3 5j9 

~Sp 



'l + r'-^-^JX{l,r\0^ 



2r 



l-r^ + ^-J\il,r^O 



— Sp 



l+r^-e-VA(l,r2,0 



K' 



n'^\il,r^,0 



for {m,n) = (-1, 0), (0, -1) 



K' 



n^mlJX{l,r^O 



for (m,n) = (-1,1), (1,-1) 
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